


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1985-03 


Analysis of incompressible cascade flows 
using state-of-the-art computer programs. 


Genskow, Paul E. 


http://ndl.handle.net/10945/21148 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


mi whaeaAny 
NAVAL £21 {SADUATE SCHOOL 
b ONTEREY, CALIFORNIA 93343 

















NAVAL POSTGRADUATE SCHOOL 


Monterey, Galifornia 





THESIS 


ANALYSIS OF INCOMPRESSIBLE CASCADE FLOWS USING 
WSTATE-OF-THE-ART'' COMPUTER PROGRAMS 


by 


Paul E. Genskow 


March 1985 


Thesis Advisor: M. F. PLATZER 





Approved for Public Release; Distribution is Unlimited 


T220185 








SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 
READ INSTRUCTIONS 


REPORT DOCUMENTATION PAGE BEFORE COMPLETING FORM 


2. GOVT ACCESSION NO.| 3. RECIPIENT'S CATALOG NUMBER 





1. REPORT NUMBER 














5. TYPE OF REPORT & PERIOD COVEREO 
Master's Thesis 


March 1985 


6. PERFORMING ORG. REPORT NUMBER 


8. CONTRACT OR GRANT NUMBER(2) 






4. TITLE (and Subtitie) 
Analysis of Incompressible Cascade Flows using 
"State-of-the-Art'' Computer Programs 













7. AUTHOR(s) 


Paul E. Genskow 


4 


9. PERFORMING ORGANIZATION NAME AND ADDRESS 10. PROGRAM ELEMENT. PROJECT, TASK 
AREA & WORK UNIT NUMBERS 


Naval Postgraduate School 
Monterey, California 93943 





t2. REPORT OATE 


March 1985 


13. NUMBER OF PAGES 
64 


15. SECURITY CLASS. (of thle report) 


UNCLASSIFIED 





ti. CONTROLLING OFFICE NAME ANDO AOORESS 






Naval Postgraduate School 
Monterey, California 93943 


4 MONITORING AGENCY NAME & ADORESS (If different from Controfling Office) 












1Se. OECLASSIFICATION/ DOWNGRADING 
SCHEOULE 


16. DISTRIBUTION STATEMENT (of thia Report) 


Approved for Public Release, Distribution Unlimited 


DISTRIBUTION STATEMENT (of the abetract entered in Block 20, If different from Report) 





18. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reverse side if neceesary and identify by block number) 


Incompressible Cascade Flows, Computer Programs Analysis 


20. ABSTRACT (Continue on reverse elde if neceeeary and identify by block number) 

Two ''state-of-the-art'' computer programs, Q3DFLOW and a standard Cebeci 
boundary layer program as applied to aerodynamic flows in a cascade are 
examined. Q3DFLOW is a quasi three-dimensional finite element turbo- 
machinery program. Only the incompressible blade-to-blade module of this 
program is used, The Cebeci boundary layer program ts for incompressible, 
two dimensional, and constant viscosity flows. It is used to calculate 


DD , Rather 1473 = EDITION OF ! NOV 65 1S OBSOLETE 


SN 0102- LF- O14- URITY CLASSIFICATION OF THIS PAGE (When Dela En 
ue Mag fo ] SECURITY CLASSIFICATION OF THIS PAGE (When Dota Entered) 





SECURITY CLASSIFICATION OF THIS PAGE (When Data Entered) 


boundary layers on a few isolated airfoils and airfoils/blades in a cascade. 
Experimental data and other analytical computations are compared for both 
programs. Additionally, the theory and operation of the Cebeci boundary 
layer program are presented. 


S’N 0102- LF- 014-6601 


2 


SSS ap a pt al 
SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) 


a 


Approved for public release; distribution is unlimited. 


Analysis of Incompressible Cascade Flows using 
State-of-the-Art Computer Programs 


by 


Paul E. Genskow . 
Captain, United States Marine Corps 
B.S., lowa State University, 197 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 
from the 


NAVAL POSTGRADUATE SCHOOL 
arch 1985 


ay LIBRARY 
YDLEY KNOX fete ecHOOL 


STGRADU 
NANT EREY ; CALIFORNIA 93943 


ABSTRACT 


Two "state-of-the-art" computer programs, Q3DFLOW and a stan- 
dard Cebeci boundary layer program as applied to aerodynamic flows in 
a cascade are examined. Q3DFLOW is a quasi three-dimensional finite 
element turbomachinery program. Only the incompressible blade-to- 
blade module of this program is used. The Cebeci boundary layer 
program is for incompressible, two dimensional, and constant viscosity 
flows. It is used to calculate boundary layers on a few isolated airfoils 
and airfoils/blades in a cascade. Experimental data and other analytical 
computations are compared for both programs. Additionally, the theory 
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l. INTRODUCTION/BACKGROUND 


Two important parameters of interest in the field of aerodynamics 
of airfoils are lift and drag. Lift is defined as the force generated by 
an airfoil normal to the mean free stream direction of flow. Lift can be 
approximated by integrating the pressure distribution around an airfoil. 
Drag is the component of force parallel to the mean free stream direc- 
tion of flow. Drag on airfoils is usually defined to consist of two 
forms, skin friction drag and pressure drag where the total drag on 
the airfoil is usually defined as profile drag. Drag on airfoils is more 
difficult to calculate and predict than lift. 

Flow separation occurs when the fluid “particles” fail to follow the 
surface of the airfoil. Separation is sometimes confused with the 
phenomenon of stall. Stall can not occur without separation but separa- 
tion can occur without stall. Separation is usually a three dimensional 
effect and incorporating this effect into lift and drag calculations on 
airfoils which are two dimensional is a very difficult task. Prediction of 
drag on airfoils on a purely theoretical basis is difficult even when the 
flow has not separated since transition of the flow from laminar to 
turbulent significantly affects the drag. 

Laminar flow is distinguished from turbulent flow visually by the 
macroscopic motion of the fluid. Laminar flow appears smooth to the 
eye while turbulent motion is very erratic to the eye. The physical 
mechanisms that induce this motion are different. Generally speaking, 
most flows of interest start laminar (steady) and transition to turbu- 
lent (unsteady). The start and finish point of transition is very diffi- 
cult to predict. This transition point or range is usually empirically 
determined. 

The unsteady nonlinear three dimensional Navier Stokes equations 
define the macroscopic motion of a Newtonian fluid. A Newtonian fluid 
is a fluid in which the stresses are related to the rates of strain of the 


fluid particle and the coefficient of viscosity. Air is assumed to be 


governed by these equations. However, until recently there was no 
oractical method of solution to these equations. Due to this difficulty 
the Navier Stokes equations were reduced to two general sets of 
equations: 

1. Inviscid/potential 

2. Viscous/boundary layer. 

The inviscid equations are further simplified by assuming steady 
flow, no external body forces,and no energy transfer to or from the 
fluid. With the additional assumption of uniform approach flow the flow 
is irrotational and hence permits the introduction of a velocity poten- 
tial. The inviscid/potential flow solution is of limited value because it 
predicts no viscous drag (D'Alembert's paradox). However,in combina- 
tion with the Kutta-Joukowski condition the inviscid solution predicts 
the lift of airfoils with sharp trailing edges at small angles of attacks 
very well. The condition originally proposed by Kutta in 1902 and 
independently by Joukowski in 1906 [Ref. 1: pp. 390-395] asserts only 
that the upper and lower surface flows meet smoothly at the trailing 
edge. This condition is sometimes referred to as a zero load condition at 
the trailing edge, meaning there is a pressure equalization at the 
trailing edge. 

The Kutta condition is a practical solution to the problem that the 
inviscid flow solution is singular when the flow negotiates a sharp 
corner (that is the velocity goes to infinity) for an arbitrary circula- 
tion. Placement of a rear stagnation point eliminates this problem and 
introduces a unique solution to the inviscid equations. Experimental 
evidence supports the theory of circulation and the Kutta condition on 
airfoils with relatively sharp trailing edges at small angles of attack. 
However, in certain flow regimes the inviscid/potential flow solution 
with the Kutta condition fails to adequately explain real world flow 
phenomena, such as flows around wings at high angles of attack and 
practical turbomachinery blades in cascades. Here the effects of 
viscosity can no longer be ignored. 

With the effects of viscosity included, the unsteady Navier Stokes 


equations are difficult to solve. Simplification is accomplished by 
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discarding higher order terms. These simplified Navier Stokes equa- 
tions are now termed boundary layer equations. The boundary layer 
concept , first introduced by Prandtl in 1904, defines a very thin 
region next to the surface of an object where all significant velocity 
gradients normal to the streamwise direction occur and hence viscous 
forces are generated. Outside this region all viscous’ forces are 
assumed negligible and the inviscid/potential solution is assumed to 
adequately model the flow. The boundary layer equations can predict 
the drag on airfoils since most drag is due to skin friction caused by 
the viscous forces in attached flow. The problem is made more difficult 
by the need to compute both laminar and turbulent boundary layers and 
to determine the transition point from laminar and turbulent flow. The 
reader is also reminded that the boundary layer concept is applicable 
only at sufficiently high Reynolds numbers. For a more detailed 
discussion of the boundary layer concept see Schlichting [Ref. 2] and 
Cebeci and Bradshaw [Ref. 3]. 

Turbulent flows are difficult to predict and model. Alleviation of 
the unsteady nature of turbulent flows is generally accomplished by 
expressing the pertinent variables of the unsteady three dimensional 
Navier Stokes equations as a time averaged quantity and a fluctuating 
quantity and then time averaging the equations term by term. These 
time averaged Navier Stokes equations produce nine additional unknown 
entities that can be correlated to a stress tensor. These extra turbu- 
lent unknowns are termed Reynolds _ stresses. However, all nine 
Reynolds stresses are not considered significant. In fact, only one 
Reynolds stress is usually retained in a two dimensional analysis, which 
is the case for the Cebeci program used in this thesis. This Reynolds 
stress is the shear stress that is in the normal plane to the surface in 
the streamwise direction. It is considered the turbulent counterpart to 
the laminar shear stress. This particular Reynolds stress is sometimes 
modeled as the product of an empirical “eddy viscosity’ and the local 
velocity gradient. 

Computational fluid dynamics is becoming a more important field as 


memory capacities grow bigger, computers get quicker, and numerical 
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algorithms get better. Accurate prediction of flows minimizes the 
amount of wind tunnel testing and helps the designer to make better 
designs. Two ‘state-of-the-art’ programs for predicting flows and 
boundary layers in a cascade the Q3DFLOW-81 and a CEBCAS are exam- 
ined. The Q3DFLOW-81 is a turbomachinery program developed by Dr. 
Hirsch and associates of the Free University of Brussels [Ref. 4] and 
[Ref. 5]. The blade-to-blade portion is used to calculate inviscid 
incompressible flows in a cascade. Inviscid results are presented for 
several airfoils/blades and configurations. Experimental data, analytical 
solutions or other computational flow data are compared where appli- 
cable. The CEBCAS is a standard Cebeci two-dimensional, incompres- 
sible and constant viscosity boundary layer program modified to give 
graphical output for boundary layer and skin friction distribution and 
calculate losses and drag coefficients for airfoils/blades in cascade. 
The CEBPLT is also a standard Cebeci two-dimensional, incompressible 
and constant viscosity boundary layer program modified to give graph- 
ical output for boundary layer and skin friction distribution and calcu- 
late drag coefficients for isolated airfoils. CEBPLT and CEBCAS are 
used to calculate boundary layers on several isolated airfoils and 
airfoil/blades in a cascade. Experimental data is correlated when avail- 
able and reasonable to do so. Additionally, the theory and operation of 
the Cebeci boundary layer programs (Chapter III. and Appendix A 


respectively) are presented. 
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Il. INVISCID FLOW COMPUTATIONS IN A CASCADE 

Physical flow is three dimensional and unsteady, especially in an 
axial flow compressor. Simplification of this complex flow in a 
compressor is usually approached by § splitting it into two two- 
dimensional flows: 

1. The meridional flow 

2. The cascade flow (sometimes referred to as the blade-to-blade) 
Solution to the flow in the cascade plane is input for the meridional 
"through flow’ solution, which makes the accurate prediction of flow in 
the cascade plane very important. Gostelow [Ref. 6: p. 86] has attrib- 
uted systematic advances in blade design to the ability to calculate 
flows in a cascade. 

Inviscid flow computations in a cascade are indeterminate for prac- 
tical turbomachine blades. The outlet condition, usually the outlet 
angle, is specified and at present is found experimentally. For isolated 
airfoils the Kutta condition gives meaningful results at small angles of 
attack for practical shapes (that is airfoil shapes that have relatively 
sharp trailing edges). However, most practical turbomachinery blades 
have trailing edges that are rounded so as to avoid stresses and manu- 
facturing problems. Frequently, laminar flow will separate eo the 
leading edge and reattach turbulently downstream. This further 
complicates the modelling process of the flow. Using only an inviscid 
model it is very difficult to predict the flow. Miller and Serovy 
[Ref. 7] examined five different trailing edge hypotheses to determine 
the deviation angle for cascades using the inviscid model. None of the 
five trailing edge conditions tested were considered acceptable as a 
general rule such as the Kutta condition is for sharp trailing edges. 
This is understandable since the outlet angle is influenced by viscous 
effects and on a round trailing edge the determination of a rear stagna- 
tion point is not as easy as on ae sharp ttrailing edge. Gostelow 


[Ref. 8] reported in a particular example solved by exact potential flow 


theory that a change in rear stagnation point location of as little as .3% 
of the chord would change the outlet angle by as much as 10 degrees. 
Also, the precise location of the rear stagnation point at some incidence 


is not necessarily the “correct” location for some other incidence. 


A. Q3DFLOW-81 (BLADE-TO-BLADE MODULE) 


The Q3DFLOW-81 is a turbomachinery program for calculating flows 
in an axial flow compressor or turbine. This program uses a finite 
element method. The blade-to-blade module was used to solve the 
inviscid equations for flow in a cascade. Substantial runs were made 
with the Q3DFIOW-81 program for various compressor blades (including 
even sharp trailing edge airfoils) in a cascade. In order to compare 
with experimental data each configuration was run twice, one time with 
outlet angle specified by the experimental angle and another with the 
Kutta condition option. 

The major problem encountered in running the Q3DFLOW was in 
generating an acceptable finite element mesh. The finite element mesh 
generation appeared to be very sensitive to the placement of coordinates 
used to define the shape of the airfoil/blade, especially around the 
leading edges and trailing edges. Observations of running Q3DFLOW 
(blade-to-blade) indicate a need for a close one-to-one correspondence 
between lower and upper coordinates of the surface of the airfoil, say 
within one to three percent of chord (i.e. the streamwise direction) for 
the major portion of the airfoil and closer for leading and trailing edge 
if rounded. The more curvature in the airfoil geometry the smaller the 
spacing between adjacent coordinates needed to define the geometry to 
ensure that flow computations are for the right shape. However, too 
many points to define the leading edges or rounded trailing edges will 
generate erroneous mesh elements. _A general rule of thumb is to use 
the minimum points needed for the eye/brain to generate the correct 
airfoil shape by mentally connecting coordinates. 

Four parameters ILWIN, IVPIN,ILWOUT, and IVPOUT of CARD NR 


B-03 of users guide [Ref. 9] are very useful in solving finite element 
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REPRODUCED AT GOVERNMENT EXPENSE 


mesh problems. These parameter modified the mesh around the leading 
and trailing edges by forming a triangle element instead of the normal 
rectangular shaped elements(Figure 2.1). Two other parameters that 
are sometimes useful are NSTUP of CARD NR B-O1 (number of stations 
in the upstream) for solving leading edge mesh problems and NSTDWN 
of CARD NR B-01 (number of stations in the downstream) for trailing 
edge problems (Figure 2.1). A ‘last ditch maneuver’ is to input 


different coordinates to define the airfoil/blade shape. 
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Illi. THEORY OF CEBEC! BOUNDARY LAYER PROGRAM 


The governing equations (Eqs. 1, 2 and 3) of the Cebeci two 
dimensional boundary layer program are obtained from an engineering 
simplification of the two dimensional Navier Stokes equations by 
neglecting small terms due to the known empirical fact that the 
displacement thickness is very small compared to the characteristic 
length of the airfoil. These equations are then assumed to be appli- 
cable to turbulent flows by including only the Reynolds shear stress 
and neglecting the two Reynolds normal stresses that are obtained by 
time averaging the two dimensional Navier Stokes equations. Reynolds 
stresses are essentially empirically obtained and can be expressed in 


many different forms, depending on the theory used. 
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The eddy viscosity formulation is used in the Cebeci program to 
model the Reynolds shear stress. This concept was first introduced by 
Boussinesq in 1877 [Ref. 2: p. 578] in which the turbulent shear 
stresses are assumed to be similar in nature to the laminar shear 
stresses and are therefore proportional to the velocity gradient (Eq. 
4). Cebeci and Smith [Ref. 10: pp. 255-257] have formulated a two- 
range empirical eddy viscosity model. In the inner region the formula 
Is similar to the Prandtl mixing length formula modified by a Van Driest 
damping parameter and intermittency factor that accounts for the tran- 
sitional region. The outer region is modelled as proportional to an 
empirical constant,the velocity at the edge, the displacement thickness, 


and also to the intermittency factor. 
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REPRODUCED AT GOVERNMENT EXPENSE 


Assuming constant density and viscosity and using definitions as 
stated in Eqs. 5 and 6, Eq. 1 becomes Eq. 7. For laminar flows b is 
equal to one. Also b may be redefined so that Eq. 7 is useful for 
axisymmetric laminar flows (Mangler Transformation) with small tran- 


sverse curvature effect [Ref. 3: p. 48-50]. 


3 5 E@case 
Beale. Eq. 6 
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Since the rate of change of local static pressure in the y direction 
is simplified to equal zero (Eq. 3) and the velocity at the edge of the 
boundary layer (that is when y is equal to the displacement thickness) 
iS a function of x only, then the Bernoulli equation (Eq. 8) may be 
used to determine the rate of change of local static pressure in the x 
direction, as the value in the free stream and surface are considered 
equal. Differentiating Eq. 8 with respect to x gives Eq. 9. (Note that 
Eq. 9 is now a differential vice partial since the local static pressure is 
a function of x only.) The local static pressure is now specified as a 
pressure coefficient, velocity ratio, etc. from the _ inviscid/potential 
solution as a boundary condition for the boundary layer equations. 
With Eq. 9 ,the simplified momentum equation in the x direction (Eq. 7) 
and the continuity equation (Eq. 2) reduce to a system of two partial 


differential equations with two unknowns u and v for laminar flows. 
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The Falkner-Skan transformation (Eq. 10) defines a dimensionless 
coordinate n and Eq. 11 defines a stream function » which consist of a 
dimensionless stream function f in x and n. A necessary and sufficient 


condition for the existence of a stream function is that continuity (Eq. 
2) be satisfied and hence Eqs. 12 and 13 must be satisfied. This 


i 


REPRODUCED AT GOVERNMENT EXPENS 


transformation along with the definition of a stream function and use of 
the chain rule for differentiating reduces Eqs. / and 2 into one partial 
differential equation with one unknown stream function in yn and x for 
laminar flows with no known analytical solution . For similar flows this 
transformation changes a system of two partial differential equations 
into one ordinary differential equation with an unknown stream function 


as a function of the independent variable n. 


n= ge y Eq. 10 
v=(uyx)?F (x, 9) Eq. ll 
un $Y Een le 
va~ Qu ag, 1S 


Similar flows are laminar flows that satisfy Eq. 14. An example of 
a similar flow is the flow over the flat plate. In this example the 
pressure-gradient parameter, m is equal to zero. The definition of m 
(Eq. 15) may be derived by differentiating Eq. 14 by x ,substituting 
Eq. 14 back into this differentiated equation, and then solving for m. 
The concept of similarity is not as important now as it was before the 
advent of the computer. It allowed for practical numerical solutions of 
certain flows even though very few of the flows can model any real 
world phenomenon. For further understanding of the concept of simi- 


larity see References 2 and 3. 
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max Sue cas 


Derivation of Eq. 22 is a simple matter of putting u , v and their 
partial derivatives in terms of the dimensionless stream function f by 
using Eqs. 10-13 and the chain rule. (For example Eq. 16 gives Eq. 
17, where prime denotes partial differential with respect to n). Eqs. 
17-21 show the results of this. Substituting Eqs. 17-21 into Eq. 7, 
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REPRODUCED AT GOVERRIMENT EXPENSE 


rearranging , expanding, and multiplying the equation by x and 
dividing by u,squared, and then substituting in the pressure gradient 
parameter (Eq. 15), Eq. 22 may be obtained for non-similar laminar 
flows (b=1). Now b is a function of n and must be written as Eq. 23 


to account for turbulent flows. 
un Fe- (aux) GESR Eq. 16 
usu Bees 17 
va $ (ex) + (PY?) - Cun)? 9h Eq. 18 
Qu. dus psy, OF Eq. 19 
gu (sss ya a 20 
gry Se. oer 2 


pmsl pp smL1- (FF) -x(r'96-#"96) Beja 22 
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This third order nonlinear partial differential equation (Eq. 23) is 
written as a system of three first order partial differential equations 
(Eq. 24-26) by defining two new variables script u and script v func- 
tions in n (Eq. 24 and 25 respectively). Script u is a nondimensional 
x-component velocity but script v has no direct correspondence to the 
y-component velocity. The three unknowns of this third order system 
are script u, script v, and f. Recalling Eq. 1/7 and the fact that the 
velocity is zero at a solid surface of an airfoil (n=0), the stream func- 
tion f and its partial derivative with respect to n is equal to zero. 
Hence two of three boundary conditions needed are Eq. 2/ and 28. 
The third boundary condition (Eq. 29) comes from the velocity at the 
edge (n= displacement thickness) being modelled equal to the velocity at 


n equal to infinity for a given x location. 


an 
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— Eq. 24 
one Eq. 25 
(bu) 54+ m[1-&?)]-x (UQe- BF) Eq. 26 
F(x,0)=0 io. 27 
ux ,0)=0 Eq. 28 
ux ,n)-0 Eq. 29 


Cebeci uses central finite differences of truncation error of second 
order in n on the net rectangle shown in Figure 3.1 and averages at 
the midpoints of the net. The central finite difference of the first 
derivative formula (Eq. 32) used in Cebeci's program can easily be 
derived by use of a Taylor series expansion of a point x+h and a 
point x-h (Eq. 30 and 31 respectively,see Figure 3.2). Eq. 32 is 
obtained by subtracting Eq.31 from Eq.30 and then solving for the first 


derivative. 
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Figure 3.1 Boundary Layer Grid 
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F(xth)=f (x)+hf (thi gher order terms Eq. 30 


f(x-h)=F &)-hF (x )+higher order terms East 


p(x) Pbth) fob), 94,2) Eq. 32 


The O(h*) in Eq. 32 represents the sum of the remaining terms of order 
h? in the series. Eq. 29 gives an exact value of the derivative if all 
terms in the infinite series are retained or if h goes to zero. Since 
Eq. 29 is truncated by discarding all terms of order h?* this finite 
difference formula is of truncation error of order h*. The spacing 
between grid points (h) is usually relatively small and the error 
incurred by this approximation is usually insignificant for most engi- 
neering purposes. Central finite differences of the first derivatives are 
applied at the mid points of the net. Twice the hx and hy distance are 
defined as kn and hj (that is the x distance between i-1] and i grid 
location and the n distance between j-1 and j grid location respec- 
tively,see Figure 3.1). Rewriting Eq. 32 tn terms of subscript 
sequence grid location numbers i and j of the spanwise and normal 
direction respectively at the midpoints and accounting for the fact that 
all derivatives are partial gives Eq. 33 and Eq. 34 respectively. For 
example, Eq. 33 is the partial derivative of f in the spanwise direction 
for any constant normal grid location (j). Using Eq. 33 to approximate 
Eq. 24 and 25 averaging script u and script v at the midpoint in the } 
direction , and rearranging, Eqs. 35 and 36 are obtained. In a similar 
manner and with use of the central difference approximation for the 
partial derivative in the spanwise direction Eq. 26 its transformed from a 
partial differential equation to an algebraic equation and will not be 
presented here due to its lengthy derivation and additional definitions 
introduced that are unnecessary but convenient [Ref. 3: pp. 216 and 
217]. 
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Figure 3.2 Points Defined for Taylor Series Expansion 


If all values of f,script u, and script v are known at some initial 
downstream spanwise grid location, say i-1 for all values of j then Eq. 
35 , 36 and Eq. (7.2.10c) of Reference 3 are a nonlinear system (that 
is the system is nonlinear since the last equation is nonlinear) of 3J+3 
equations with 3J+*3 unknowns where j=0,1,2,3...,J. Cebeci uses 
Newton's method to solve this system of nonlinear equations. Cebeci 
departs slightly from the strict application of Newton's method in 
treating the b terms of Eq. (7.2.10c) of Reference 3 to avoid unneces- 
sary complication [Ref. 11]. The solution to this now linearized set of 
equations is accomplished by a very efficient block elimination method 
which is possible due to the unique form that the system of equation 
can be put in. This method is fairly elaborate and will not be 
presented here. 

The basic idea of Newton's method can be understood from a little 
geometry. Suppose we wish to solve for the roots of this nth order 
polynomial in the independent variable x where n is greater than 1 and 


hence nonlinear. Say this polynomial f(x) is defined as in Figure 3.3 . 
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Now suppose we guess a solution to one of the roots of this polynomial 
say x1. From geometry, the tangent of 6 is f(x1)/x1-x2 which is also 
the slope at x1 that is f'(x1). Solving for x2 and writing in terms of 
sequence iteration numbers for successive approximations we obtain Eq. 
or: An exact value to this root can not be given with only a finite 
number of iterations so an acceptable convergence limit is usually 
defined between successive iterations. Cebeci considers the wall shear 
parameter that is f' at n equal to zero (also Known as script v) as the 
greatest source of errors in boundary layer equations. The conver- 
gence criterion is .00001 between successive iterations of the wall shear 
parameter (f''(0)) for laminar flows. For turbulent flows the conver- 
gence criterion is expressed as a percentage. There are many problems 
that can be encountered by this technique. For instance if the initial 
guess is "poor" this technique may converge to a root not sought or 
may not even converge at all. However, due to the nature of the 
boundary layer equations and the manner in which Cebeci applies it, 
Newton's method is “unconditionally stable" For a more detailed discus- 


sion of the theory of Cebeci boundary layer program see Reference 3. 





Figure 3.3. Newton's Method 
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IV. RESULTS OF O3DFELOW INVISCIDSBEABESIO Bes DE 
COMPUTATIONS 





Several blade shapes and configurations were run using the blade 
to blade portion of the Q3DFLOW-81 program. Several different graph- 
ical outputs are available from Q3DFLOW, but only those that could be 
compared with other available information are presented here. Results 
are presented in graphical output usually as pressure coefficients. 

A no camber,sharp trailing edge airfoil, the NACA 65-010 (see 
Appendix B for shape and coordinates used to input Q3DFLOW) was run 
at an inlet B = 30°, o = 1.5, and several as. Figure 4.1 shows the 
results of the Q3DFLOW-81 blade to blade (solid line) plotted as a pres- 
sure coefficient vs. chord. The circle and triangle points depict the 
experimental pressure coefficients obtained from NACA TN 3916 
[Ref. 12]. Figure 4.1 shows the results specifying the experimental 
outlet 8 as given in NACA TN 3916. Figure 4.2 shows results using 
the Kutta condition. 

A slightly cambered, rounded trailing edge blade, the NACA 65-(4A, 
lg5, )10 (see Appendix C for shape and coordinates used to input 
Q3DFLOW) was run at an inlet 8 =30°, o =1.0 and several a's. Figure 
4.3 shows the results of Q3DFLOW-81 blade to blade (solid line) as 
pressure coefficient vs. chord when specifying the experimental outlet 8 
as given in NACA TN 3817 [Ref. 13]. Figure 4.4 shows the results 
when using the Kutta condition. Experimental pressure coefficients 
(circles and triangles) obtained from NACA TN 3817 are also plotted. 

A slightly cambered, sharp trailing edge airfoil, the NACA 8410 
(see Appendix D for shape and coordinates used to input Q3DFLOW) 
was run at an inlet B = 44.4° and o = 1.33. Figure 4.5 shows the 
results of Q3DFLOW-81(solid line) using the Kutta condition, as velocity 
ratios referenced to the approach velocity vs. chord. The velocity 
ratios, as calculated by the Schlichting-Method and given’ by 
AGARD-AG-220 [Ref. 14: pp. 317-319], are referenced to the approach 
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Figure 4. Ms Sie LOW Results for NACA 65-010 
(Outlet 8 Given from Experiment) 
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Figure 4.2 Q3DFLOW Results for NACA 65-010 
(Kutta Condition) 
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velocity vice the mean velocity as originally given and are presented as 
circles and triangles in Figure 4.5 . 

A cusped blade (see Appendix E for shape and coordinates used 
to input Q3DFLOW) from Gostelow [Ref. 6: pp. 109,110,123,124] was 
run at an inlet B = 37.5° and o =1.01. Figure 4.6 shows results as 
pressure coefficient vs. per cent chord using the Kutta condition. The 
circles and triangles on Figure 4.6 show the results as analytically 
determined by Gostelow. 

An advanced controlled-diffusion fan blade (see Appendix F for 
shape) was run at an inlet § = 68°, % = 58.25° ,o = 1.34 and an axial 
velocity/density ratio (AVDR)=1.443. Figure 4.7 shows results. of 
Q3DFLOW-81 as pressure coefficient vs. chord (solid line) using the 
experimental outlet angle as measured in the cascade wind tunnel at the 
Turbopropulsion Laboratory of the Naval Postgraduate School. Figure 
4.8 shows results using the Kutta condition. Experimental pressure 
coefficients as measured in the cascade wind tunnel at the 
Turbopropulsion Laboratory of the Naval Postgraduate School are 
plotted as circles and triangles. 

Table | shows outlet B as calculated from Q3DFLOW using the Kutta 
condition option for all blades run. The Kutta condition option of 
Q3DFLOW uses an initial outlet B which was given as the experimental 
outlet 8 or other analytical outlet 8. Also shown is outlet 8&8 given from 


experiment or another available method. 
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Figure 4.5 Q3DFLOW Results for NACA 8410 
(Kutta Condition) 
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Figure 4.6 Q3DFLOW Results for Cusped Blade 
(Kutta Condition) 
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Figure 4.7 Q3DFLOW Results for Fan Blade 
(Outlet 8 2 Given from Experiment) 
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Figure 4.8 3DFLOW Results for Fan Blade 
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V, RESULTS OF GCEBEGC) BOUNDARY LAYER SPROGRAM 


Several isolated airfoils were run on CEBPLT. Several airfoils and 
turbomachinery blades were also run using pressure or velocity distri- 
butions from inviscid flow in a cascade (Q3DFLOW-81 blade to blade) as 
input for CEBCAS. Results from CEBPLT and CEBCAS are presented 
as graphical output displaying displacement, momentum and _ skin friction 
distributions. The velocity distribution input is also displayed. 
Transition was computed by Michel's method as given by Cebeci and 
Smith [Ref. 10: pp.332-333] or the point of laminar separation if this 
occurs first. Drag coefficients on isolated airfoils were calculated using 
a Squire. Young formula [Ref. 2: pp. 764-765]. Losses were calculated 
for airfoils/blades in a cascade as given by Roudebush and Lieblein 
[Ref. 15: p. 24]. The drag coefficient for cascades was calculated by 
multiplying the losses by the cosine 8 at infinity (i.e. the 8 associated 
with the mean free velocity) and dividing by oa. 

NACA 0012 was run using velocity ratios obtained from an inviscid 
program of the Douglas Aircraft Company at a = 0°. NACA 0009 and 
NACA O006 were run using inviscid velocity ratios from Abbott and von 
Doenhoff [Ref. 16] also at a = 0°. NACA 4412 was run using the 
experimental pressure coefficients from Pinkerton [Ref. 4] again at a = 
O° . Figures 5.1-5.5 show the results of these airfoils as displacement 
thickness, momentum thickness, and_ skin friction distributions vs. 
x/c, fraction of chord. The profile drag coefficient ,C,, and skin 
friction drag coefficient, C, ,referenced to approach velocity and 
velocity at the edge,u,are for one surface only and are given on the 
skin friction graph. The total profile drag for symmetrical airfoils is 
twice the C, value given on skin friction graph at a = 0°. Table II 
shows the comparison between calculated profile drag and measured 
profile drag as given in [Ref. 16]. 

The NACA 8410 airfoil was run on CEBCAS using the velocity 
distribution as given by the Schlichting method in [Ref. 14: pp. 


Jf 


rmEPROOUWCED AT GOVERNMENT EXPE WN@ 





NACH OOl2 AUR=O Re -oevenS 
THICKNESS DISTRUBUTION 


x10" 


a 
™~ 
wm 
ea) 
oe) 
a 
awe 
= 
a. 
an 





Ave 





Figure 5.1 NACA 0012 Boundary Layer Results 
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Figure 5.2 NACA 0009 Boundary Layer Results 
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Figure 5.3 NACA 0006 Boundary Layer Results 
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NACA 4412 Boundary Layer Results (Suction Side) 
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Figure 5.5 NACA 4412 Boundary Layer Results (Pressure Side) 
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317-318]. after referencing to the approach velocity as input at a 
Reynolds number of 500,000. Figures 5.6 and 5.7 show results as 
thickness and skin friction distributions with transition specified as 
close to the leading edge as possible without incurring a program error. 
Also shown are the velocity distribution inputs. The displacement and 
the momentum thickness as calculated with the help of Truckenbrodt s 
method for a fully turbulent boundary layer as given in [Ref. 14] are 
plotted as squares and circles, respectively, on the thickness distribu- 
tion graphs. 

The NACA 65-010 and the NACA 65-(4A, lg, 10 were run on 
CEBCAS using the inviscid pressure distribution from Q3DFLOW as 
input at Reynolds number of 245,000 and 444,000, respectively. Figure 
9.8 and 5.9 show results as thickness and skin friction distributions for 
the NACA 65-010 at an a = 6°. The drag coefficient for the cascade 
referencing to approach velocity as calculated by CEBCAS for the 
NACA 65-010 at an a = 6° was .0138 as compared to .013 as measured 


in NACA TN 3916. Table I!l shows where separation occurs according 
to CEBCAS in x/c. A "NO" indicates that separation did not occur 
over forward 95% of chord. Lieblein's equation used for losses does 


not account for separated flows or blades with round trailing edges and 


for that reason no other comparisons were made with experimental data. 
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Figure 5.7. NACA 8410 Boundary Layer Results (Pressure Side) 
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Figure 5.8 NACA 65-010 Boundary Layer Results (Suction Side) 
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Figure 5.9 NACA 65-010 Boundary Layer Results (Pressure Side) 
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TABEESH 


Comparison of Calculated and Measured Drag Coefficients 


a REYNOLDS EXPERIMENT MEAL CULATED 
AIRFOIL 
NACA 0012 O° 6 MILLION . 00600 .00562 
NACA 0009 0° 6 MILLION .00S%3 .00614 
NACA 0006 OF 6 MILLION .00510 .00610 
NACA 4412 Q° 3 MILLION .00584 .00693 


TABLE II] 
Separation Points from CEBCAS 


NACA 65-010 (Inlet B=30° o=1.5) 


Approximate x/c 


ao Upper Lower 
-3° al ¥O2Z5 
23 mes, 0235 
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Vie DISCUSSION/CONCLUSIONS 


AS me ISCUSSION 


The experimental pressure coefficients as found by the program 
Q3DFLOW agreed well with experimental data when using the Kutta 
condition or inputting the experimental outlet 8B with the exception of 
the fan blade. Q3DFLOW also gave good results when compared with 
the Schlichting and Gostelow analysis. Disagreement between experi- 
ment and Q3DFLOW for the fan blade case may have been due to the 
fact that the incidence angle run was very close to complete stall (i.e. 
wind tunnel stall) as evidenced by subsequent runs in the cascade 
wind tunnel at the slightly higher incidences. The axial velocity/ 
density ratio as measured in the wind tunnel was very high in this case 
and indicates lack of two-dimensionality. Q3DFLOW accounts for this as 
a blockage factor through the. cascade but is only an linear approxima- 
tion. These two factors may have contributed to the disagreement 
between the Q3DFLOW results and experimental data as given for the 
fan blade. Q3DFLOW may give unrealistic pressure coefficients also 
One encountered was as high as -34 . The coordinates used to define 
the leading and trailing edges may have been the cause since a judicial 
choice of the coordinates generally alleviated the unrealistic pressure 
coefficients. 

The discontinuous pressure coefficient distribution at the trailing 
edge on NACA 65-(4A,1g,)10 (Figures 4.3 and 4.4) is probably due to 
the discontinuous geometry at the trailing edge that was input to 
Q3DFLOW as a "last ditch maneuver’ to generate an acceptable finite 
mesh. Initial attempts to generate an acceptable finite element mesh 
using coordinates where the trailing edge was_ sufficiently defined were 
unsuccessful. (Figure 6.1, coordinates were obtained by blending in the 
forward 95 % of blade with a 1 % chord radius trailing edge located as 
given in NACA TN 3817.) NACA TN 3817 indicates that the NACA 
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65-(4A, 15, )10 may have possibly separated on the concave surface 
(pressure side) at as of 3.4° and 5.4° for inlet 8 = 30° and o = 1. 
Figures 4.3 and 4.4 did not show any significant difference between 


measured and calculated pressure coefficients. 





Figure 6.1 NACA 65-(4A,15,)10 Shape 
ann Round Trailing Edge Defined) 


Experimental data on boundary layers in cascades is sparse at 
best which makes it difficult to evaluate the merits of the Cebeci 
boundary layer program. A further difficulty is that most cases run on 
the Cebeci boundary layer. program using Q3DFLOW inviscid incompres- 
sible pressure distribution as input indicated separation. High pres- 
Sure gradients and low Reynold numbers on the leading edge may make 
the boundary layer equations not valid. However, the calculated drag 
coefficient for the one case that had a sharp trailing edge and had 
attached flow on both sides, the NACA 65-010 at a = 6°, did agree 


well with the measured value. 


B. CONCLUSION 


Results of the program Q3DFLOW indicate that reasonably correct 
pressure distributions may be calculated when inputting the leading 


edge and trailing edge coordinates of a practical compressor blade with 
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care. The usefulness of the Cebeci program in calculating boundary 
layers in cascades still remains a question. Use of these two programs 
in conjunction with each other will probably not give accurate aerody- 
Namic parameters such as lift and drag on practical turbomachinery 
blades in general due to separation being calculated by the Cebeci 


boundary layer program. 


31 


APPENDIX A. 
OPERATION OF CEBEC! BOUNDARY LAYER PROGRAM 


The boundary layer program as sent from Cebeci in September, 1984 
calculates total skin friction drag, total profile drag using Squire- 
Young formula, and transition based on an empirical formula given by 
Cebeci and Smith [Ref. 10] in addition to those parameters presented in 
Cebeci and Bradshaw [Ref. 3]. Additionally, the inputs have been 
modified. The current version of this program used for this thesis has 
been slightly modified again to give graphical output and to handle 
numerous data files at once. Another slight modification can calculate 


losses for cascades using an equation 44 of Reference 15. 


A. INPUTS 


All inputs are formatted by fortran format. A card will refer to an 
input line that is a maximum of 80 columns. Each card will be defined 


by columns,format, variable and explanatory comment. For example: 
col 1-5 15 variable explanatory commentary 


col 1-5 indicates the field assigned to hold 


the value of the variable 


\e indicates the format used 

variable indicates a fortran name for the 
variable 

explanatory comment gives information on the variables 


Integer format statements are right justified that is in the preceding 
example if you want the variable to have a integer value of 1 you would 


put 1 in column 5. If you put the 1 in column 3 the value would be 


o2 


read as the integer 100. The other two fortran format statements used 
in this program are F10.0 and 20A4. The F10.0 refers to a floating- 
point variable. A real number must be specified within the ten position 
field using a decimal point to separate the decimal point from the 
mantissa. Any amount of numbers right of the decimal point may be 
specified as input as long as it fits in the field. The 20A4 refers to a 
Hollerith variable. 80 positions are available for 80 characters. This 
format is used to read in the title for the output and the disspla 
routine (graphics). The last two character of the title must be $' so 
the disspla routine may read your title. 80 characters are available to 
specify the title but twenty to thirty are only recommended for the 
current character height used in the program. This may easily be 
modified for those familiar with DISSPLA ( a= graphics’ software 
package). 


t. inputs for CEBRLT 


The following cards give the format to the modified Cebeci program 
(CEBPLT) that gives graphical output but does not calculate losses for 


cascades. 
SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 
Run card Format(15) 

eK RRR IK 

col 1-5 I5 N the number of runs (data files) 
SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS$ 
Title card Format(20A4) 

KEKKKEEKEKE 

col 1-80 20A4 Title tabular and graphics title; 


last two characters must be $ and 
immediately next to last character of 
actual title, recommend only twenty 
to thirty characters for current 
character height setting 
SHSSSSSSSSSSSSSSSSSHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS$$ 


General card 1 Format(315) 
KKKKKEEEEK 
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CO ll Ome 


col 6-10 I|5 


eo} | 2 [ss 


NXT total x/c or y/c coordinates; 
number of x/c coordinates must equal 
number of y/c coordinates (i.e. 
they are paired) 

ITFS transition flag setter; 
set equal to 1 transition calculated 
by Michel's method which is given 
by Eq.(6.1.1) in Ref. 2 set equal 
to Q, transition is input XCTR 
(General Card 2 col 41-50) 

BC boundary conditions format; 
Set equal to 1 when inputting 
velocity ratios, UE, on General card 
2, col 31-40 , set to any integer 
value other than 1 when inputting 
the Coefficient of Pressure ,Cp, 
instead of UE 


SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS$SS$SS$S$SS$$$s 


General card 2 
Rk EK 


col 1-10 F10.0 


col M=ZO0s bree 


col 21-30 F10.0 


Col 31-40 F10.0 
Col 41-50 F10.0 


Format(6F10.0) 


el, Pressure gradient parameter; 
value of m at stagnation point, set 
equal to 1.0 

DET At initial normal step size,h; 
(i.e. is in the neta direction), 0.01 
is a good value 

Gi Variable grid parameter, k; 
defines the ratio of successive step 
sizes in the neta direction , 1.14 ts 
a good value,the higher the value 
the fewer calculations performed in 
the neta direction 

RC Reynolds number of chord. 

XCTR transition input; 
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x/c coordinate of transition from 
laminar to turbulent, required if 
ITFS is 0 (General card 1,Col 6-10) 


SSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 


Daca card 2 
KKK KR RK 


col 1-10 F10.0 XC 


col 11-20 F10.0 YC 


eo 1-30 7F 10n0 UE 


2. inpiits sion CEBCAS 


Format(3F10.0) 


x/c coordinate; 
starts at stagnation point, or leading 
edge,must be in increasing size 

y/c coordinate; 
must correspond to the x/c coordi- 
nate 

Velocity ratio; 
that is the velocity at edge of 
boundary layer divided by velocity 
at infinity, must correspond with 
location (x/c,y/c) and with BC 
assigned a value 1 (General card 1, 
col 11-15), Coefficient of Pressure 
may be inputted instead if BC is a 


value other than 1 


Inputs to CEBCAS are essentially the same as CEBPLT except for the 


run card. The run card is modified as follows to input other parame- 


ters needed to calculate losses and drag coefficients for cascades. 
SSSSSSSSSSSSSSSSSSSSSSSH$SHSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSSS 


Run card 

tok ink ek 

col 1-10 I5 N 

col 11-20 F10.0 BTA 
col 21-30 F10.0 BTA2 
col 31-40 F10°0 SIGMA 


Format(110,3(F10.0) ) 


the number of runs (data files) 
the inlet B (degrees) 
the outlet B (degrees) 


e) 


Do 


B. OPERATION OF CEBPLT AND CEBCAS 


First you must obtain the appropriate fortran program. Link to the 
AERO disk and send appropriate program to your disk. After you got 
a copy of CEBPLT or CEBCAS on your disk you must compile the 
fortran program(s) on an appropriate fortran compiler. FILEDEF'’s (file 
definitions) must be defined for data input and for tabular output. An 
input data file must be created as defined. Execute DISSPLA at a 
TEK618 graphics terminal for graphics output. For specific details see 
appropriate computer manuals at the Winston Churchill Computer 


Center. 
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APPENDIX B 
NACA 65-010 (SHAPE AND COORDINATES) 
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)10 (SHAPE AND COORDINATES) 


APPENDIX C 
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APPENDIX D 
NACA 8410 (SHAPE AND COORDINATES) 
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CUSPED BLADE (SHAPE AND COORDINATES) 





4h 4h db tb tp 2h 


th tt Sh db ib db tb SE Sh ob 4b db ob dE db de gh 2b de ce at at at 
+ 4 


tt 44} 

+} | a § g 4h <> 
st 4b 4? BY 4? 
aT 2 SY g 4} st TH=NNIFNODDOMNIRODRNONR§ 4} 
th ETE OOO KOM tT eR NOR eM EMOEmMo +4 
4 $$ FUT OM Nm QO mM NIV EIM Fine op 
4} 47 TXT UDOOVOVAOOOVOCUDOVOU0OOO 4} 
+P tt UGE Pp © © ee ee eo em em we em em wm wee eo ee ow os 
4> th de H*OOOCTOODDOOOVOCUVDOOOOOCO}} 
st te As} ) ies et |e Jae a Ye | 4 
$< 4 ue 4t $> $$ 
sow Sh U St +} 3} 
s}Ws: Ode tt 4% 
te Se UN SE Sb Sb tb Sb db db db dt tb dE SE OP St SE db 4h ot db SE ah ah db ae $f 
43 Wdt th ott an 
ce > 3b X J oT g st 
st Ub 4S Us dp +¢ tt 
7% St BH SE ON DUD FT ONOCO KH AMQnNeHDaDO+4} 
se Wot U 4p $F OMOR KN NNNNMN TOWNE OTOMMO 4} 
se mm Sb tb OU 2p I= XM ONRNANOWODOTeRER NO} 
HL UNG OD DOO HRM NNMM SF IFN OKRA WKC OY 
$$ & 4} tt Mit © © © © © we eo eo ee eee ee ee 8 le ot 
{$e 2} +¢ H$ODODOOVOODODCOCOOCOCCOO0O00 4} 
S}Oss tbh fe vt 
4¢ Ur tt th $f 4} 
ft Sy es oY 2 at 
7% Ove se 4s 3 
Tt OU St th db Sb St th de 4b Sb oh dt dh Sb Sb Sh th Sb Se SP ab ct ot oP $$ 4b 4b 3t 
4} $¢ 4} oe. +b 
3b St tt st it 
se Osho 4b Sh $$ 
west fot MDMODENNBR OF TOOUNINO 4 
3 © 4} + 6 UWE MOMDANNGU MOR WOWS 
A oe Fs med dal Les | ma~NITRK OOD == OOUR UUIM of 
FUE NTE OD OOOO aa mmm tO OOOCOO} 
eM, OS OP ot ee e 3b 
4¢ mast St 4} HFOOOTOVVOOOOVUNDVUVOVOUND OOF 
4% st UL 4p 3t t 
th mein de = tf t+ 
4 (% S$ J 4} $$ +f 
4 b $$ UW) 4b % tt 
S$ Wot Sb th ib th dE de db ob db Sh dh tt dE de Se OSE OE St SE db ch Se oe ap $} 
So Rob st tt + 
St O 3 W 4h +} +} 
shor Use dt tt 
3¢ Ott Ut Sf %%$ VOOYVUYUNUYUNOMWYUIANOUNDO- 
$$ ? ef FRPUOFSAHNIMNYR FTOOMOONO4 
$$ uu) 2b te WJ eb Ce MIOIOIAFNANDON-WONKR WO: 
2OwW YUNGOOCKCOMMNANMN~ FC SINORAT L208 
se & Sh oh} KM sh 8 © © oe oe Oe gp POP eo oe ee ww oe @ 4} 
% Jt tf 6H OC OOOO O COCO OOOOOC OCHO =F} 
3h Ud Sp 4¢ 0 Sf 4} 
4¢ rf i td 4 tt 
st +t $} st +t 
2 4$ 4h t¢ 4} 
tg ob TP 4b db Sh dh 4b tb Sh SE 4b db db dt Sh OSE Sh Sh Sh db db SP ot Sh de se 4p 
| i rt 4} 
sf a g 2} st st 
2} sh be Sh SE NMPFMVOK DHOHNM FNMOKNR DCO =m 4} 
$$ St eth ont ot omnt emt omnt met emt ont mt ot OJ NJ Sh 
+} 4¢0CSE ii tt 
2 2} ST g +t $ 


t 
3h Fh th th Sh db tt de de gh th at 4h db db 4S th 4h gh ab de sh dt dE de at th ab ap 


60 


REPRODI'CED AT GOVERNMENT EXPENSe 


APPENDIX F 
FAN BLADE (SHAPE) 





61 


10. 


1) 


12 


ioe 


LIST OF REBBRENCES 
He Ad K., Principles of Ideal-Fluid Aerodynamics, Wiley,New 


Schlichting H. , Boundary-Layer Theory, McGraw-Hill Book 
Company,New York, 1968. 


Cebeci T. and Bradshaw P., Momentum Transfer in Ze ndary 


und 
Layers, Hemisphere Publishing Corporation, Washington, 197/. 


Hirsch Ch. and Warzee G.,_"A Finite Element Method for 
Bee ales Calculations in Turbomachines ,_ Transactions of 
noes ME, Journal of Fluids Engineering, v. 98, pp. 403-422, 


Hirsch Ch. and Warzee G., ‘An Integral Quasi 3-D Fin 
Element Calculation Program for Turbomachinery Flows ,ASM 
4 oun aa! of Engineering for Power, v. 101, pp. 141-148, J 


Gostelow J. P., Cascade Aerodynamics,Pergamon Press, 1984. 


Miller M. J.,Serovy G. K., "Deviation Angle Estimation for 
Axial-Flow Compressors Using Inviscid Flow Solutions , Journal 
of Engineering for Power,April,1975. 


Gostelow J. P., Trailing-Edge Flows Over Turbomachine Blades 
and the Kutta-Joukowsky Condition, ASME-Paper 75-GT-94, 1975. 


Author _unknown, Computer Program_ for Turbomachinery 
Flows/Finite Element Method/User's Guide,11 May 1983. 


Cebeci, T. and A. _M. O. Smith, naiysls of Turbulent Boundary 
Layers, Academic Press, New York, 1974. 


Cebeci T., Keller H. B., “Accurate Numerical Method for 
Boundary-Layer Flows_1!!: Two-Dimensional Turbulent Flows" , 
AIAA Journal ,Nov. 1972. 


Emery J. C., Erwine@. Ro Silene L. J., Systematic 
Two-Dimensional Cascade Tests of NACA_ 65-Series Compressor 
Blades at Low Speeds,NACA TN 3916 ,195/. 


Emery J. C., Erwin J. R., See: M., Two-Dimensional Low 
Speed Cascade Investigation of NACA Compressor Blade Sections 
peace Systematic Variation in Mean-Line Loading, NACA TN 


62 


14. 


Scholz N.,translated and revised by Klein i from a riginal 


K eo 
German Aerodynamik der Shaufelgitter, AGARD-AG-220, 1977. 
Roudebush W. H., Lieblein_S., Theoretical 42 ss Relations for 
Low Speed Two-Dimensional Cascade Flow,NACA TN 3662, 1956. 


Abbott |. H. and von Doenh 
1959. 


Off A. E. ,» PReory of @¥agg 
Sections, Dover, New York, 


Pinkerton R. , Pressure Distributions over the Mid-span Section 
of the NACA 4412 Airfoil, NACA Rept. 563, 1936. 


63 


INITIAL DISTRIBUTION LIST 


Defense Technical |Information Center 
Cameron Station. 
Alexandria, Virginia 22314 


Library, Code 0142 
Naval poet raduate School 
Monterey, California 93943 


Chairman, Dept. of Aeronautics 
Code 6/ 

Naval Postgraduate School 
Monterey, California 93943 


Professor O. Biblarz 
Code 6/Bi 

Naval Postgraduate School 
Monterey, California 93943 


Commanding Officer 

Attn: Captain P. E. Genskow 

Naval Aviation Rework Facility 
Marine Corps Air Station | 

Cherry Point, North Carolina 28533 


64 


Copies 


10 


















a ‘* . aie et 7 | ae! 
y “etre . est °. a4 . Ye 43, f ¢ : ahs . : : : | . : 
a rok To Jet £25 ay, Ff bers A ci LAY ® 
af git ise ts ! ieee, Riko sos: eine He ESoe Rar pega thesG25977 
Sb sah ee a tel RA AP A NC Eecaie ee aR NS! ita tig 
OS aera : 


; | e flows | p- 
4 ; ; 1 HI Aging bas : 
Co] } : } : : Ws f yf 1 * ad : 
DRS ee K-. : . it fi : bts ; Ache ie Y the wee 2 P a ha 
eae . Pi Seon y 4 AS What gh J + 4 aA bs I 3 


* Mm WAN a Sea 
oo 






ri 
ays 





Sins 
“eg tent 'e 


















































































































































































































































































































































































































































































, ' ' 
f . 1» ' 
i P A 1 ae 
Mantes 7%} ? : , A : 
= ie es i 
Hatt fairs AS acraN TE ar ot oy 3 2768 002 02570 i . 
AN AT Ol 2 0 Aj Ll f i EG Mb ass CU a ta dg t 
ee ig patra eae aS : ; ARY : is 
Be dnd i aCe CaM erated iat teat ge oA OX LI | ea 
KAtpn 8h Oi 0 : APA AAS BTAVUGE patna Ah pdf ASP Why “ DUDLEY KN ' ; : 
ta 4 rae F ed 88s oad Mi: aa td f | eS Ot @ ; Mp of f ; ‘ F f y 
Mts ’ be ah gi betes ad A PUB: fein mater hate 8, Rey 7 : ' ' ' 
i ¢) ay RE AG. Mah Fang \, F i Rs 4 % ware ier . é - 
BO roe EDEL Sed EE nity Ke MN a isa ; ca at ve 
f bead ua yt bk 3 th im, - U 
As ijn BR Ye sayd a 
PAS NEL Est A) tude ET vs * 
VA # irik Lie ee 4 i : “fe 
(ax Tuli | aoe , : F , 
ae 
1 Eigen ge ‘ : ' ' 
er KS ey J ' 5 7) «6 ° ‘ . hing 
9 Ro. a Lah ig os ’ KY ie 4 ‘ ci these 1 : 
: CPU far OE ay a ee Bic. ial gti , ge ‘ le Ig ‘ o o° tae 
=H Hdl, Faye a) i  frate tr, gt hat ogee P Bee: t al i pe ‘ : ‘ 
v : OS fig if t+ ‘ Sums ¢ ‘ e ‘ ¢ oa , ‘ 
Pitt ; Seen ites, rn ¢ong Mri?! yy ' ; ea Ae ; , rr 
fa os é ah onad tos og yp fee $i Pag Spe aU Amar iar A Ta F ' ea a@ 
‘ ; egy! thet at ede ate { 46 © ie at t an ‘ ot 0 ne: ve 
f - rs An i P : y * ae es A it Vat, a al fay f a ae r! M4 1 Lei 
4 ITER by . D Tike ve LS hk “a a iF foie, ee { ? 5 J ‘ ; 1 = ‘ a" 
rd A BA ou anche Ge ADA NS ‘¥ a Peat oe ‘ ' weer 
ae Loge pi rtD . ‘ Ma py a ; ‘Soue : : 
Oe af a2, Ia ee: Cote Me ‘ oi 1¢ ' 
Am celia ae shh ch ae A hile ; 
° ee fe an eae ‘ ' 
: we 1 ps4 : 11 
Ui Or @ ‘ be : ae 
pages 9.4 er 2 Say re , 
#4 ° 
er, tke 2 Oe ts 4 . I 
otrdag 8 WARD eniy hl a al 4 1 r oP. 4 
sake Shae ae a) off a, ' it foo4 
Nw aie FF onl ee, Ci J ee a ' 
MR NG toe, Bho ‘ke el Gone ety po th « ‘ 1 ' 
Odes Ob Bol td stigg ge) mca ar ‘ ' 
. CMs ab ‘ F ‘ 1 ° f 1 | i 
Fd had OL, @ Ode ee Ae) . ' ° eae ‘e4 ‘ 
aloe tier es t Site ry *! “ff \ po erie oy es ‘ 
Bye th dime? 006. Pang i fe egin te Kees “her, Aee z Ue ee at 2 ere eat 4 ‘ 
SP rt de iteloo wind spieth ss beh a or fet Oe af Kris bh aase dA ¥ a) a eae 5 ; : ! eS é ba P . ‘ ‘ fs Peet 8 es 
- ot Ky, F ot thhe ve * fia Tr : ") : ‘ ' ‘ * ° 
SG eh oat ab 9 ft chad Vener, Ag? om st THe : Fe a : : ; ; rae d 
Abdi ob oe sah . oa ae ? rs ] fete 4g * f eo. 
oP AS nip shoe @E ad encerts sg Morac ot et Ie ‘og ' Lc ela ear "teapda Hie a, tia ‘ ' v4 
Wodt thew ae bed ye Fat, “ip sy Pes ne wn to4 Re ¢y cn 4 ' e ree 
x FAL whi oles wy SEE cere te aie * lati ie, haa » if - 4 e Pt are fs is LI AEC ¢ oN 1 
PO BT oh we grap Potincene . cays OOM Seung ¥ oned jy Pe a) eS a . y Ee ee eves ‘FFs tgn, eee os . 1 os ie ' 
tab enor 5 Wd: AYP red Fibs TH ES Seb kid ee Trice a meri ems ens 1 ae Lea) Pg , ae te 
oF ass * ts dak ge oh o3tee Pd at i ‘ pete tae nat Di ns Pari Bie , ‘ e U ‘ * ry ‘ o¢ 
a Mile co Z, et ee Hrohnd Pe 7 7 ee in 2 hat J " Le ‘ 1 ee ste ‘ . ° bars 
itatclua ig Lad cet E tH OF ha cae nig ee fa Woe eres Ue Te ery es a i lieu ‘ : ' A 
nap pagiben ll EAST ED shud ; ‘ oot Pty » th ae on¥aas gh "4 earl te Soe ee ar * 5 ‘ 6 Gare 1 ’ 
SH in doe A ah lh at Ps Faugeras Ud ae ee ey. | r i dare ’ a isnt *fy twee 4 f ere te ee. ' ao » ' Ps ' } ‘ . . 
Apt hs, Mules bie 4 Se Be abe. A r : : * Dnapeee of Poteet Ce yl OF ee ier Ore a i ter fier aa 1 ‘ Ie ry ' 
Wee =. f, a 7 i ey fs RAE es gig ene r) ‘ : « 4 4 ‘ rs 
ated Oe Pre $ ese ay AU Ja Rae ln SOO a er ' Hoe , ' ' 
Olah A or iey, Sl ee TAP Bed wl af Fmt AA ? P s - Lair A @#at es ee: a Fa ee gas eos ‘ ‘ ' 1 
chdinw eat eat Piper, Shs PE hate ate FOS Wohi oan ‘ ‘ #4 we WSL tal ST iC Sri Jar ga e ta 1 en ae 1 ? » 
oe we FW bp RA tis oy) : Pe Sie JD é Hes Clit wee afar tare . die Of as 11 fee fey 
Tons mee tip. P) Ol ome shia 4 Ve “fh Fapiites of 288 4 fo f 1¢ 44 e 2 1 ‘ 1 ‘ 
ON é bate th AP pati Pat gua ey eae ¥» y . eng! a? : gy DAE tary | if hie} i) ' 1 ry . 
1G: OM Dh ie MA ® ae? ei i tree eT ae aeen A | ort ieee CURURCI |W herr ° ‘ 
F shy AOR omy os ir eS, Lae | thy AO ear 4 ra ' 
“4 ‘ Pie ON SCAR hes , Hae . ‘ ¢ 
° E eta Pt gy ape site ¢ei ire Sear ON, 
; * 6 ate oh od6 hd Cabs Ne LP aap : 
ae Sik f°ars "Akh aay wie 2h 7 i ‘ 
SF 02d ot ce #02 eb 01d oe od Ae pe RG ISAS hea ain ene t . 
A ov PU Cour we age , te ee ae boos 
ies tr ei © Fay + pe Vows? 8 E, ee er © mF bee, 
CN oat et Le awd Pete ¢ fi Ws ae fefafehaw Q 
SoS 00 i Hel bas ales tel LCR We tii Ziutsy a5 
“o ig BE PA aes ag Poca dik é F itye% of, hing ale Cpe ogee 
Wie 2806 ob oe od se tat Or eh + 
via wns phe AR dw Put ce Wd Oth & Peptrd 5 - 
ee ere * : ne MEF F, wae ¢ 
80M Mie TK d 
Westie 2) & Ss 1 eeu ma Cau 
TOW bso S r ree Dag Mak AIP IS saat rye # t6 
PE eH ot 6% Wh Matra Fy tron mre ain gb naa OS F 
Decoaate eb PyoD tan 0 Ait 16 8 0 Sel ogg poe 7, a ‘ 
Vi Dak ieee TEP Settee hia eek Le fee heh Sd vey Fee 
6 ef whee ato why oot sigs ee ee Cae hy 
fhe wid of s6@ taf ° Fh ate hg ig 
a Steen at o“ 
vw Pmiete Ned ere eins 
Lt ? Suess 6 ig 1 
4 ae 
of 09 Ss 
4! fe wee 
wL 
CV PRIA yet 
vr ; J ark, 
a\e.r # ener mann ety te : hike 
F ON Re (agai . 
A hy ete nm, ap ene, wey =a 
~ «<O™t ¢ dehy 
‘he "SNE ° 
Scar Rees y Ae Sin oe, my « ay “iv ewes, < Yroiweys ‘2 
Bl take On WEE Ris ve *% ~ * a ae 3 
ES Hee the 1 oe me Preaek Se it (2 rarer vee 
Sete 880 88 ga Qiao mesh Th we a bay e UE NIN OMe oy 
73 HE vee BMA As wy, “wa ® y eFy ‘ Ty OB Saree, i 
«ESE M, 8p Fig hen “SVG we y Wig wy der 
is) lee eee FeP pe, heey te Sy » 
Dar tee et MO) 42. 4 : SEAT N Ne Tw oy ae ‘ 
Oe Som, Mae. bat ieee Phe S nen ak ae be ora en ' 
TLI* Pay aine ye Sire Whe bine > te SL MPLA, De yr ek be tap ee 2 1 ’ 
kitten: babe. *. ah ah He tg ns Fapae yt ‘ ms eee a : is See wy A , a) - 
4 nde, te oh ba hot fara 9 Pers 9 Rein He dn oy he bel ce Lok at a * seal 3 4. Re Saag 
Tee Oey emg artyncen Me wei Sows te ree gtk i Ay ey NK on te Ate Pe Mirah Cad oti ale eens ME CORN 40 NNO Gates. 
Seen eg bey tees tht | Wy Marae orgs ay bait Pe et SO ee Wate tk Sora BAS sASis Se ety fet vee YNTy 9 +t a 
My th feaehice tt “ahh IS NEV IB Ye Qa dy ys Ae ry, 1 Sig Me my > 4 * 2a foen oy - e:. m™Je A 8 tee ns 1 VM re50 @ S a ; 
¥ em, Ae ram: Ae LM Me my amy SEE VV Bese #4, Neve te, ka en ere av “Ae aan, i LNG Meets wawrere ia 
aa town ens i bah heats he. wear, - Catt, Sn Ye SS 17 ee Ve PUN, ty, St Pe~y aie aes . uh x > OBA tee 
othe - c “itp & "e POPES yt om soy Ap rern &. ay adhe ET okt te ¥ 4 Wee o's &,* i <n Tree : on 4 h se | 
dda pote a 8 ee teh he ts MORO AE RS? i ™ Sarre ALY snipe os wre phe bs a wr Near wan A 

2 , Crthade hated th teed ree ee is et WR sRin betta ee SY ie lh bk ~EYiw op . 
rae de ch disldh sang sa Ktidindy tn Setints ants Toh ne Wale Reta Rua teams TR tee ee CL ANe on wut : ' 
dn Ne Sy tee mrt enedy Stee TS VY, Heo Wty sy . Wrmade a, Nites i i mS * ts ha Rees 5 * : 

2 were oN Ale Wee Se Ya ert % Yo in ee a : ae i 
a us ewN yy & 3 Ve heey uM, Vere rwsanas Tete oP as . ay aa, a utane * a rv 8 : : 
as > e \ Reha crt hy PAN wre Wiayt> Gitte oy. Leen koe anaes Fy acs . “A 
Pe Pay “ay = NOT chge td Sua NHN NDS * Paix VPM, whee ery “¥'S Uy ee gry ORG ' 

fi Md SER sel tats’ mars te ng eat ae bs MOPED yQus tA SUG ey. mw Oo + es Te Ne tne 
Le erate my, n, Ee pes cP i wen he os ure Sm § he OD 8 oe we Saineeg ; 
Rt Ae alae eh ets Ne: aig hs SSDS UN Oy ot NFS ae eae Her, “= e ee 14 see t atl nh . 
eet toe ete “ar Naryinn eS * aintats Sesh as re OU rey, SF By RN Qe tip ng ¥ omen, ens ¢ 
eel: Soe ese en eer weed SEAN ea hes RL ah itee en ee Nd a 
bhp Minha de Ade mda ee ALY emati re to Fy YS 6 2 tune yay ak Ms Wray Soh WK ano “HECIDRORGT ir i 
KR N Navy wouee se eps amt Face ze. Ab tetas, NEP GY Sy tat ot Sp Fes qitrya 5 ‘* t Ve enwAn 
WP EAT hp a te iptlenta len te Lt AO eve Be ty: Se UEt WIA Ne Situ . 4 we Th eis 
1 Fim. * ge BO? ew fey rey: TS Oe dry My Bp Mam, wy, & &* eae ? 

Fite hy 4 cer eran bitin Sea ta es , 1 whey & te i : 
7S Deg ou | bale ia de an TY te iy lwrng oe WH ae be WI UDET) 816 von Lp ty ns ' = 
eine he Tee ie eM *. VFEN AVE Te FH omny “ve, & 19 re Wie 2d TY ; ‘ ery aD aes ‘ 

adel chelate So ey my: NTP Tq WW 9 & cHrdofend Then ee ey ety : 
Sturt hal bedapeeh Crhe Saat can % ° TE eID ie aves ie ¥ I A alee oh A Mk % 
phe bial - VA ely Swarr enatas 7 + bol 7 5 ey ep wae es a { t + ; & As ee © *¥ hig a i a ay 
oe Loeeestre Weta rag nted Ee : nee: Malt =< “Wok Woe BE Ay Pate Ria SP a ek Fa ot, ES Ke ts re Vee ase ee Syne . 
ate yg Ee ee Avene meh ae? e ele ars 21S eS. : OS Se het : 
ae My area tat tat ee gah ZVI Lary $ 
nla aia he Te FPP ey he 2G ey: Yiere COHN “e ' . 4 
Shek th Ae SS WI, tahy a4 Oy ’ 

Me ta ae LO. OS hee re Uma ” “WMP aesy ¢ 
bat Nes’ AVES A NG THT I k Nie G40 tage ne ere © koe h hs ‘ 
mar ee ayy 3 PE e729 Fee 8 ht a Bee ee UU A ony wa £ Pauses ‘ : 
PN My Pep y ean a ot EEG De ye rasan 6.4 ‘aye p"® SUV Pes N a uey oe ngs Gok Win toy : 
wy Garem, ary Ph hpi ler te doh hy “Athol det hb Pt) tee Tee a PEW UERA te +. . 
bart la nag et ee SIEM Foray ALE try, bys 4) IES Pak we os Se ‘ 7 epeee 
FPS he Bh gr tig tats lan TAY MO Hs oy ADSL De wpa Spb beak Ce i 3 : 
Mahe ne LS a ’ NEN d Keeani i ae 3 \e > : bee a 4 
Sed Fa ae Pay. eae ne AY Y BOer ear Ee iy PREF Sing SFR ORB FAL 6) le ’ iy ety i d any if » ee eee atace: 
let ee TT ae ‘ NIU WeR GS Eyer cork) OPe TR RET Aang Ws ed DO Ree, SAS “~hes . eed | o SFU DONE he wy 
RS See ~~ ¥ ital ata 2h ety 971 10 Fy Ong ty ON Were Bog, ; < : ‘ “tte elegy 
P-beeiin Ma he te Pott Lh is Meo ky OE REY tarpey 0s 4 VAY Wel, wae WYER ob ey @ kre say, a Cire : 
5, veyron: Ss eet ite Tost } ny atlbehd hated P89 Wrdae age 
oe Wty a, my O18 42 « SONA abe epi 94 iny SR DIere: rE Vey #9 Leh tase oe ow tet F 
Tete ore SUES salyrn wy Ny C4 Madan hg LAAT La] hy Taare ae : 
Reality Sat PY SOs Way tpn BP REN Ome 8 ay hy © Ay! SN W"Hy WU Ea ae ty po ; 
Hid, Sate in CEES ry VA he 1 ey Ay ey SHAN TE LEY S600 6 Se iy 2 yi PPR AR Ey fay ote 8 sony Ser, 
Cee AVI Sy Ae Cartes 29 m0 vty bladed beta Wot tie os Wego ere eRe ey FIRED Mire eo : 
Cures BES ANS my & Pager ay he, ING Ni rte Bp me eye ns LO WE WaitD nc whe 
*& pm. pede anlar oy WOK ervey eae tts ny» 
og Fe ay PAT Bet SON TUD psy CVG ey Renova gta ton ed 3 
Reg eS. PEK Bion EP WRONG V9: SOG WY Weer y ctynd yeh DL, Oh RP SO ds de ot A We MES wre any 
eres baad et Wb Ae PP AE ema one a) PUN Pe Py oe ye eRe Me les Aa nabee ber pity °t ‘late Pitan et 's 
OCF ENTE Fh, tn ge 9 My ey REE TM, ww Kat Sw UEUPIEE A TEL Nee ara ia E 
bie A ak iE Mi tte tata ete aia Na Mee” Bae. 
BEVRAL Brian wy: a NI 8 ot Mr ty Me ete nee His eek Crs 
eNO PER 84-0 1@ eG D8 baer Sey is ARE ban SY 
% HR JAAR RAs yn 
Py t 7 


A ee 
be On ee bf 
Ly A hehe et 
\O Wis 


Me 8 Oy 
My te ae O°) freon oe oe, ‘ 4 
© 8 Berea crm *Usie.g 
Mach de 2 2 ry 
VREROE ISI we: vr 








UNG § oF © wre 

eet ANG ath Mer won Ste 

PRAY He & 01064, ow i tris, Phy YUVA mo > mw ; 

reste SLT een Pai a Rete aaa vi rick vs ree Ny aye a? we 
tarhey tates TAY ROS PRONE yma Le MM HOLD Le 
* Faye ube » ib wd 20 5 cane = 

Fe ee RODD, eg 


